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Zn  this  paper  me  describe  some  results  ob- 
tained over  the  past  few  years  in  the  areas  of 
estimation  and  control.  The  point  of  this  research 
is  to  find  subclasses  ->f  systems  that  possess  suf- 
ficient structure  to  allow  us  to  determine  optimal 
solutions  or  computationally  feasible  algorithms. 
Zn  particular,  we  utilise  several  results  from 
group  theory  (both  for  Lie  and  finite  groups)  to 
obtain  our  results. 


between  these  structural  extremes.  Zt  is  our 
feeling  that  results  of  this  type  not  only  pro- 
vide useful  solutions  to  certian  problems  but 
also  shed  considerable  light  on  the  implications 
of  structure.  That  is,  we  can  begin  to  answer 
questions  such  as  which  results  really  require 
linearity  and  which  can  we  extend  to  systems  with 
somewhat  less  algebraic  structure. 


ZZ.  Finite-Dimensional 


The  results  we  review  in  this  section  were 
obtained  in  collaboration  with  James  Lo  and 
Steven  Marcus.  For  thorough  treatments  of  these 
and  related  results  and  for  further  references, 
we  refer  the  reader  to  [1-8] . Consider  a dynami 
cal  system  described  by  the  following  equations! 


Although  in  principle  many  control  and  estima- 
tion problems  have  been  solved  in  extremely  general 
settings,  the  solutions  are  often  unacceptable 
because 

(a)  The  "solution"  does  not  actually  provide 
a design  (as  in  the  so-called  Kushner- 
Stratonovich  solution  of  the  nonlinear 
filtering  problem) . 

(b)  The  resulting  design  is  much  too  complex 
to  be  evaluated  for  a specific  system 
(as  in  solving  extremely  high-order) 
Riccati  equations) . 

or  (c)  The  design  is  computationally  too  complex 
for  on-line  implementation  (as  with  many 
nonlinear  filtering  algorithms  or  the 
control  of  large  interconnected  systems) . 

Often  these  difficulties  arise  because  of  the  ex- 
treme generaltiy  of  the  problem  being  solved.  The 
approach  of  the  research  described  in  this  paper 
has  been  to  restrict  attention  to  subclasses  of 
systems  and  problems  in  order  to  exploit  the  alge- 
braic structure  that  they  have  in  common.  Our 
motivation  has  been  that  there  must  be  a middle 
ground  between  "linear"  and  "nonlinear"  or  between 
"small-scale"  and  "large  scale”  — i.e.  that  while 
in  both  cases  the  former  categories  may  be  too 
restrictive  for  many  applications,  the  latter  are 
far  too  large  to  provide  many  useful  results.  In 
this  paper  we  describe  several  results— all  using 
group-theoretic  concepts — for  classes  of  systems 


(given) 


d£(t)  - F?(t)dt  ♦ Gdw(t) 


where  £ e Rn , * t ^ or  R*3* , and  w is  a standard 
Brcwnian  motion  process.  Such  a system  is  often 
called  a bilinear  system,  and  systems  of  this  type 
arise  in  a wide  variety  of  applications  (see 
[1-3, 5-7]).  One  very  important  application  arises 
when  x(t)  is  a 3x3  direction  cosine  matrix,  repre- 
senting the  orientation  of  a rigid  body  with 
respect  to  inertial  space  (see  [6,7]  for  details). 
Here  the  A^  are  all  skew-symmetric  and  the  £i 
represent  components  of  the  angular  velocity  of 
the  rigid  body  (in  this  case,  of  course,  an  equa- 
tion of  the  form  (2.2)  holds  only  for  spherically- 
sysmetric  bodies;  otherwise  the  equations  contain 
nonlinearities,  reflecting  differences  in  the 
body's  moments  of  inertia). 

The  estimation  problem  we  wish  to  consider  is 
as  follows i we  observe 

dz(t)  -H£(t)dt+  h1/2dv(t)  (2.3) 


where  R > 0 and  v is  a standard  Brownian  motion 
process  independent  of  w.  Given  s(s) , 0 < s ^ t 
we  want  to  devise  a recursive  procedure  for  esti 
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HtlBf  X(t)  — i.S.  V*  wish  to  cospitt 

x(t|t)  - stx(t)  |*(a) , 0<a<t]  (2.4) 

Problems  of  this  typa  arise  in  rigid  body  attitude 
estimation  problems , when  we  have  a strapdown 
system  that  provides  measurements  of  angular  ve- 
locity (see  [6])  or  when  we  have  an  inertial  plat- 
form (see  [6,7]). 

In  order  to  understand  this  problem  let  us 
consider  the  simplest  case  — N-n-k-1.  In  this 


r a t Ai  f 5U)d»  i 

-|e°  . ° J: 


The  solution  then  is  clean  let 

p(t)  - ft  C(s)ds  (2.6) 

0 

and  design  a Kalman  filter  to  compute  the  optimal 

estimate  (£(t|t),  p(t|t))  of  the  augmented  state 
(C  (t)  ,p  (t) ) . Then 

A_t  Kfi  (t)  | 

x(t|t)  - e 0 B(e  1 |r  (s) , 0<s<t)x0  (2.7) 

and,  since  p is  Gaussian  when  conditioned  on  x, 

we  can  compute  x(t|t)  as  a function  of  p(t|t)  and 
its  covariance 

^ A_t  + A p (t|t)  + i A20(t|t) 

x(t|t)  - e 0 1 21  (2.8) 

G(t|t)  - E((p(t)-p(t|t))2]  (2.9) 

where  O can  be  precomputed. 

The  key  question  is  how  this  simple  result 
extends  in  the  multivariable  case.  Note  that  if 
all  of  the  A^  commute,  we  can  write 


x(t)  “ exp 


| Apt  + J C1(S)ds|xQ  (2.) 


and  in  this  case  one  can  readily  extend  the  scalar 
result  (see  [1]).  When  the  A^  don't  comnute  the 

situation  is  far  more  complex,  and,  in  fact,  in 
order  to  obtain  results,  one  must  determine  the 
manner  in  which  the  A^  don't  conthute.  This  leads 

to  the  introduction  of  some  notions  from  the 
theory  of  Lie  groups  and  Lie  algebras. 

Definitions (1)  Let  L be  a subspace  of  kxk  matri- 
ces (over  C) . lisa  Lie  algebra  if 


(A,B)  « AB-BA  t L 


va,beL 


(2.11) 


(2)  Let  S be  a set  of  kxk  matrices.  Then 
there  exists  a smallest  Lie  algebra  {s}^  that 

contains  S.  This  is  called  the  algebra  generated 


(3)  Let  1 be  a Lie  algebra.  The  Lie  group 
G associated  with  L is  the  smallest  group  (under 
matrix  multiplication)  containing  exp(L)  VI£  L. 

(4)  Let  1 be  a Lie  algebra.  A subalgebra 
S is  an  ideal  if  - 


IL,S]  - { [L,S]  |lc  L,  seS>  C S 


(2.12) 


If  S is  a set  in  L,  we  let  {s}^  ^ denote  the  ideal 
in  L generated  by  8. 

(5)  A Ida  algebra  1 is  nilpotent  if  the 
lower  central  series  of  ideals 


i°-l 


Ln+1  - [L,L“l 


(2.13) 


terminates  in  {<>}  for  some  n.  Rote  that  L is 

abelian  if  and  only  if  L1  * {o}. 

We  refer  the  reader  to  [1,5,6]  and  in  parti- 
cular to  [9,10]  for  physical  and  mathematical 
motivation  for  the  introduction  of  Lie-theoretic 
concepts  into  the  problem  of  analysis  of  bilinear 
systems.  We  only  note  the  followings  we  asso- 
ciate two  Lie  algebras  with  the  system  (2.1) 


L - (A^,...^ 
L0  " 


(2.14) 


Then  if  GQ  is  the  Lie  group  associated  with  LQ, 


x(t)  c eAt  Gq  xQ 


(2.15) 


For  more  detailed  controllability  results,  we 
refer  the  reader  to  [9,10]. 

We  have  already  stated  that  if  L is  abelian, 
we  can  obtain  a finite-dimensional  optimal  esti- 
mator for  x(t) , and  in  fact  that  estimator  con- 
sists of  a Kalman  filter  followed  by  an  exponen- 
tiation. More  generally,  we  obtain  the  result 
[6,8] . 

Theorem  2. 1»  Suppose  LQ  is  nilpotent.  Then 

x(t|t)  can  be  computed  via  a finite  set  of  (in 
gegeral  nonlinear)  stochastic  differential  equa- 
tions (driven  by  the  innovations 

dv(t)  - dz(t)  - nf(t|t)). 

The  basis  for  this  result  is  as  follows t if 
Lq  is  nilpotent,  then  one  can  readily  show  [6,8] 

that  every  element  of  x(t)  is  a finite  sum  of 
terms  of  the  form 
-t 

J e*5(s)ds 


♦ (t)  - e n(t) 

where  n(t)  is  a finite  Vol terra  series 


(2.16) 


(k^f  • • • 


n(t)  • 


£ /...  / E 'J 
3-1  0 0 * 

■v  yv-"V<,vn-«) 


where  the  kernels  are  separable  - l.e.  the  w’s 
are  of  the  form 

• Yq (t) Tfj (0^1  • • • "T j (®j)  (2.18) 

Theorem  2.1  then  follows  from  the  more  general. 

Theorem  2.2  ([6,8])t  The  estimate  t(t|t)  can  be 

computed  with  a finite-dimensional  set  of  sto- 
chastic differential  equations  if  n is  given  by 
(2.17)  with  all  separable  kernels.  This  result 
includes  far  more  system  than  those  described  by 
(2.1).  We  refer  the  reader  to  (6,8,11)  for  exam- 
ples and  more  on  systems  with  such  Vol terra  series 
representations. 

The  proof  of  tills  result  is  quite  lengthy 
and  complicated,  and  we  refer  the  reader  to  (6,8) 
for  details.  We  indicate  only  a few  of  the  cri- 
tical ideas.  The  separability  condition  is  es- 
sential to  have  finite-dimensional  realizability. 
Consider  one  of  the  simplest  cases 


♦ (t)  - / w(t,s)£(s)ds 
0 


(2.19) 


Using  a conditional  version  of  Fub ini's  theorem, 
we  have 


4(t|t)  - J w(t,s)£(s|t)ds 


(2.20) 


At  first  glance,  it  appears  that  we  have  to 
smooth  the  entire  £ trajectory  in  order  to  estimate 
4.  However,  if  w is  separable,  then  we  can 
realize  (2.19)  by  a finite-dimensional  linear  sys- 
tems which  when  combined  with  (2.2)  leads  to  a 
finite-dimensional  Kalman  filter  that  can  compute 
(2.20).  The  situation  for  higher  order  Vol terra 
terms  is  more  complicated  but  similar. 

Now  define  the  conditional  smoothed  covariance 

*ij<erl»<Vt)  “ El«i(a1)4i(o1|t))(5j(o2)-!j(o2|t)) 


A critical  property  of  9 is  that  it  is  separable. 
TO  see  the  importance  of  this  consider 


t a a 

4(0-  / //  cx  (Oj)  C2  «V  C3  (o3)do3do2da1 


(2.22) 


In  this  expression,  write 


W “ t^i(ai>“^i(aJt)1  * (2.23) 

I 

One  of  the  terms  that  results  from  this  is 


t o.  a 


4(t)  - Iff  tci«V"V0ilfc>1 


0 0 0 


iW-v°2|t)  )t3(o,!t)d03d03d01 


Then  it  is  easily  shown  that 
\ai  12 


4(t|tJ-  Iff  pl2<<Wt,*3  (a3l  t>dff3d02da1 


0 0s 


-/JL  Pl2(ai'flr2't,P3(°2lt,d°2dCTl  (2.3 


where 


P3«..  - /? 3< 


Again,  it  is  easily  seen  that  the  separability 
of  P32  allows  us  to  augment  £ with  a finite  set 

of  variables  in  order  to  determine  a Kalman  filter 
that  computes  $(t|t). 

One  of  the  other  terms  in  (2.22),  (2.23)  is 
t ax  o2 

M(t)-  f f f ^1(o1|t)^2(02|t)t3(a3|t)d03da2d<^ 

0 •'o  0 


Taking  the  differential  we  obtain 

i o. 


du(t)  - £,(t|t)X(t)dt  ♦ f f f 
1 0 •'o  J0 

d (?x  (a3 1 1)  t2  «t2  1 1)  C3  (03 1 1)  Jda3do2doi 


(2.26) 


(2.21)  X(t)  •//  tjtaJtK^oJtldO^  (2.27) 


Using  Ito's  differential  rule  on  the  second  term 
in  (2.26),  computing  dX(t)  in  much  the  same  man- 
ner, and  iterating  these  arguments,  we  obtain  a 
finite-dimensional  set  of  nonlinear  differential 
equations  (driven  by  the  innovations)  for  the  com- 
putation of  y (see  (6) ) . The  general  result 
follows  from  arguments  along  these  lines  and  from 
related  properties  of  conditional  Gaussian 
processes. 


Thus,  ve  see  that  this  result  combines  soma 
relatively  simple  properties  of  linear  filters 
with  the  algebraic  structure  of  nilpotent  bilinear 
systems  in  order  to  provide  one  of  the  largest 
classes  of  finite  dimensional  optimal  nonlinear 
estimators  that  in  some  sense  can  be  viewed  as  a 
"natural"  extension  of  the  Xalman  filter  [1,6] . 

3.  Fast  Estimation  Algorithms  for  Certain 
Harkov  Chains 

Let  xtk)  be  a finite  state  Markov  process  with 
state  set  S * and  with  probability 

vector  p(k)  (here  p(k>A  ■ ProbtxlkJ-S^)) . The 

probability  transition  equation  for  the  process  is 

p(k+l)  - Ap(k)  (3*1) 

where  A is  an  KxN  stochastic  matrix.  Suppose  we 
take  a sequence  {y(k)}  of  noisy  measurements 
of  the  process.  These  measurements  are  charac- 
terised by  the  conditional  distribution  ( if  y is 
discrete)  or  density  (if  y is  continuous) 


y(k)  • (xtk)  ♦ v(k))l  mod  ■ (3.6) 

where  v is  a sequence  of  independent  ^ -variable* 

independent  of  w.  Let  n(k)  be  the  distribution 
for  -v(k)  (here  "minus"  is  Interpreted  modulo  N ) . 
In  this  case,  several  straightforward  calculations 
yield  the  filtering  algorithm 

*-l 

p(k+l|k)A  - g CjPOtlk)^  (3.7) 

( , n0c)l-y(k)p(k^k~1)l 

p(k|k)  A (3.6) 

where  all  indices  are  to  be  interpreted  mod  H. 

Let  us  identify  Z^  with 

where  s is  an  indeterminate  and  we  multiply  with 

■ 0 

the  rule  s ■ 1-z  (hence  this  is  a group  isomor- 
phic to  Zjj).  Then  any  function  ft  Z^  + C is 

represented  in  the  form 


C(V.i>  - P(y(k)-v|x(k)-Si) 


we  also  assume  that,  conditioned  on  x(k),  y(k)  is 
independent  of  all  other  x(j). 

The  problem  of  computing  the  conditional  dis- 
tributions p(k|k)  and  p(k+l|k)  for  x(k)  and 
x(k+l)  conditioned  on  y(l) ,...,y(k)  is  concep- 
tually trivial,  and  the  solution  is  given  by 

Diffusion  Update 

p(k+l{k)  ■ Ap(k|k)  (3.3) 

Measurement  Update 


p(k|k)A  - C(y(k),i)p(k|k-l)i 


IV 

£c(y(k)j)p(k|k-l)  . 
j-1  3 


Although  the  derivation  of  (3.3) , (3.4)  is  a tri- 
vial application  of  Bayes*  rule,  the  implementa- 
tion of  these  equations  can  be  computationally 
nontrivial.  For  example,  in  general  we  must  per- 
form N2  multiplications  in  carrying  out  the  dif- 
fusion update.  The  question  we  wish  to  consider 
in  this  section  is  the  investigation  of  special 
classes  of  finite-state  Markov  process  for  which 
the  estimation  procedure  (3.3),  (3.4)  can  be  per- 
formed efficiently.  The  results  wo  discuss  in 
this  section  and  other  related  topics  are  treated 
in  much  more  detail  in  (12-15]. 

To  motivate  the  approach,  suppose  x(k)  C Z^ > 
the  group  of  integers  {0,1,...,N-1;  with 
addition  defined  modulo  M.  Furthermore,  suppose 

x(k+l)  » (x(k)  + w(k)J  mod  N (3.5) 

where  w is  a sequence  of  identically  distributed, 
independent,  random  variables  with  values  in 

and  with  distribution  £.  We  also  assume 


F - £ f.X1  (3.9) 

i-0  1 

We  equip  the  set  of  all  such  functions  with  the 
usual  definitions  of  pointwise  addition  and  scalar 
multiplication  and  with  a convolution  product 


**  t 

- Ev 

t-0 


where 


- £ 


(3.10) 


(3.11) 


This  set,  denote  by'C[Z ] is  called  the  (complex) 
group  algebra  of  , and,  with  this  notation 

(3.7)  becosms 


p(k+l|k)  - C«p(k|k) 


If  we  also  equip  ClZ^I  with  the  pointwise 
product 


(l9)i  « Vi 

and  the  evaluation  map 


(3.12) 


(3.13) 


.<o  - £ 


(3.14) 


then  (3.8)  becomes 


<3.15) 


Motivated  by  these  observations,  we  consider 
a store  general  problem  (for  a far  more  general 
case,  see  [14-15]).  to  do  this  we  need  the  fol- 
lowing ((16)) 


XOt)  • n*y(k) 

Y(k|k)  - A(k)p(k|k-1) 


(here  we  regard  y(k)  as  an  element  of  C[Z  ]). 

If 

examining  these  equations  we  see  that  the  diffusion 
update  (3.12)  consists  of  a cyclic  convolution, 
while  the  measurement  update  consists  of  a permu- 
tation (3.15),  a pointwise  multiplication  (3.16), 
and  a normalization  (3.17).  Mote  that  if  we  wrote 
(3.12)  in  the  vector  form  (3.3),  A would  be  a cir- 
culant  matrix. 

The  preaence  of  a cyclic  convolution  in  these 
equationssuggests  the  possible  use  of  transforms. 

Tor  f 6 C[Z  1,  define 


We  also  endow  C [X]  with  the  pointwise  product  and 
the  evaluation  map 


j2ir/n 


Note  the  inverse  transform  relation 


Applying  transforms  to  the  algorithm  (3.12) 
<3.15)-(3.17)  yields 


C^(P  (k+l|k))  -NC1(C)Ci(p(k|k)) 


C.tt(k))  - y 


(4)  Two  representations  and  T^  are  equi- 
valent if  there  exists  an  invertible  matrix  S such 
that 


H-l 

C.(Y(k|k))  - £ C. (A(k))C.  , 


Ci(Y(k|k)) 

NCn(Y(klk)) 


Examining  these  equations  we  note  two  things s 

(1)  Comparing  (3.12),  (3.16)  with  (3.21), 
(3.23)  we  see  a striking  "duality". 

In  the  CTZg)  domain  the  diffusion  up- 
date is  a N convolution  and  the  central 
step  of  the  measurement  update  is  a 
pointwise  product.  On  the  other 
hand,  in  the  transform  domain  the 
situation  is  exactly  reversed. 

(2)  We  can  exploit  this  duality  to  re- 
duce the  overall  computational 
complexity  from  0(N2)  to  0(n  log  n) 
with  the  aid  of  the  fast  Fourier 
transform  (FFT). 


(3.33) 


Consider  a random  process  xOc)  and  a sequence  of 
observations  y(k)  on  X described  by 

x (k+1)  - w(k)x(k)  (3.34) 

y(k)  - v(k)x(k)  (3.35) 

where  w and  v are  independent  sequences  of  inde- 
pendent random  variables  on  X.  Let  K be  the  dis- 
tribution of  w,  and  let  n be  the  distribution  of 

▼ \ In  the  group  algebra  domain  we  obtain  the 
following  filtering  algorithm 


Interconnections 


The  results  described  in  this  section  are 
preliminary  in  nature  and  a more  detailed  treat- 
ment will  be  given  in  the  forthcoming  report  (18). 
He  also  refer  the  reader  to  results  in  (19],  (20] 
that  are  along  these  lines.  The  basic  idea  of 
this  work  is  to  consider  control  and  estimation 
for  systems  that  are  made  up  of  many  interconnec- 
ted subsystems.  The  major  problems  that  one  is 
interested  in  are  either 


p (k+l|k)  - C«p(k|k) 

(3.36) 

A(k)  - n*y(k) 

(3.37) 

y(k|k)  - A(k)p (k|k-l) 

(3.38) 

p (k|k)  - Ii£]£2 — 
s(Y(k|k)) 

(3.39) 

If  we  consider  the  transform  domain  version  of  the 
algorithm  we  find 

C*(P  (k+l|k))  - Ci(C)Ci(p(k|k)) 

(3.40) 

(1)  Algorithmic  — i.e.  the  systems  are  so 

complex  that  the  usual  design  and  analysis 
algorithms  (computer  solutions  of  Lyapu- 
nov or  Riccati  equations,  etc.)  cannot 
handle  them.  One  wants  to  devise  special 
algorithms  for  analysis  and  synthesis  that 
. utilize  the  system  structure. 

or  (2)  Implementation  — i.e.  often  the  usual 
system  designs  cannot  be  implemented  in 
real  time,  and  one  wants  to  utilise  sys- 
tem structure  to  speed  up  implementation 
(for  example  by  taking  advantage  of 
sparseness)  or  to  determine  faster,  per- 
haps decentralised,  designs. 


(3.41) 

(3.42) 


and  we  can  make  the  following  remarks t 

(1)  The  convolution  (3.36)  of  the  diffusion 
update  is  transformed  into  a pointwise 
multiplication  (3.40)  of  the  transform 
matrices. 

(2)  The  transform  version  of  the  pointwise 
multiplication  (3.38)  is  in  general 
quite  complex  (see  [15]).  However,  in 
certain  cases  — e.g.  when  X is  a 
metacyclic  group  (161  — it  can  be 
interpreted  as  a convolution  of  the 
transforms,  yielding  the  "duality"  en- 
countered earlier  in  the  case. 

(3)  In  order  to  take  advantage  of  pointwise 
multiplication  versus  a convolution,  we 
need  to  be  able  to  perform  the  trans- 
forms (3.32),  (3.33)  quickly.  A general 
procedure  for  this  is  not  yet  known, 
but  in  certain  cases  (such  as  the  meta- 
cyclic case)  we  can  obtain  generalisa- 
tions of  the  FFT. 

He  note  that via  Myhill-type  constructions,  general 
finite  state  Markov  processes  can  be  viewed  as 
evolving  on  semigroups  (12-14),  and  thus  the 
issue  of  the  complexity  of  semigroup  algebra  mul- 
tiplication is  important  in  understanding  the  com- 
plexity of  the  filtering  problem.  Work  in  this 
area  is  continuing  (17] . 
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As  in  the  preceding  sections,  our  approach 
consists  of  examining  systems  that  possess  certain 
structural  properties  that  we  can  exploit. 

Consider  the  dynamical  system 

x(k+l)  - Ax(k)  ♦ Bu(k)  + w(k)  (4.1) 

y(k)  ■ Cx(k)  + v(k)  (4.2) 


i 


and  suppose  we  want  to  design  a feedback  controller 
that  produces  u(k)  based  on  y(j) , j <k  and  that 
minimizes 


j . E 


I £ lx«  (k)Qx(k)  + u*  (k)ttu(k)] 
(k«0 


(4.3) 


The  solution  to  this  problem  is  well-known  and  it 
consists  of  a Kalman  filter  to  estimate  x(k)  fol- 
lowed by  a calculation  of  u(k).  The  design  invol- 
ves the  solution  of  2 Riccati  equations  of  dimen- 
sion n-dim(x). If  n is  large,  this  is  an  extremely 
difficult  problem.  In  addition,  the  on-line  im- 
plementation also  requires  a large  number  of  com- 
putations per  stage. 

Suppose,  however,  that  the  overall  system 
can  be  decoupled  into  smaller  subsystems  by  change 
of  basis  on  the  state,  input,  and  output  spaces. 
Then,  we  can  solve  lower-dimensional  problems  and 
things  become  simpler.  The  difficulties  are  the 
following 

(1)  In  general  It  is  not  easy  to  find  such 
basis  transformations. 

(2)  In  general  one  can't  expect  to  decouple 
everything  (including  the  cost  J). 

(3)  for  on-line  implementation  one  would 
have  to  perform  the  input  and  output 
transformations  in  real  time.  In  ge- 
neral this  is  a difficult  computational 


task 


Consider,  however,  the  following  situation: 
our  system  is  made  up  of  N identical  subsystems 
numbered  from  0 through  N-l,  and  they  are  intercon- 
nected in  a circularly- symmetric  fashion.  That  is, 
x,  uy,w,  and  v are  composed  of  subvectors  for 
each  subsystem,  and  if  we  use  a subscript  to  denote 
the  subsystem  number,  we  have 


For  this  problem  we  have  reduced  the  algo- 
rithmic complexity  (by  obtaining  lower  order  pro- 
blems to  be  solved)  end  the  complexity  of  imple- 
mentation. That  is,  we  implement  the  system  as 
follows:  (1)  obtain  measurements  y^(k) » (2)  Use 

FFT  to  compute  y(k,l)i  (3)  process  this  data 
N"*l 

through  —j—  decoupled  Kalman  filters  and  op- 
timal gains  to  obtain  u(k,t)i  and  (4)  use  the 
inverse  FFT  to  obtain  u (k) . 

The  above  analysis  barely  scratches  the  sur- 
face of  what  can  be  done  with  such  systems. 

Other  issues  that  we  have  considered  and  are  stu- 
dying presently  are: 

(1)  Fast  algorithms  for  stability  tests, 
Lyapunov  and  Riccati  equations,  pole 
placement,  etc. 

(2)  Extensions  to  other  types  of  group 
symmetries  (with  the  use  of  the  cor- 
responding transforms  as  described  in 
the  preceding  section).  Also  the  pos- 
sible (simpler)  use  of  number-theoretic 
transforms. 

(3)  Decentralized  control  (viewed  as  a type 
of  "subsystem  low-pass  filter”). 

(4)  Utility  of  such  control  laws  when  the 
system  isn't  quite  of  this  form. 

Results  in  these  and  related  areas  will  be  repor- 
ted in  (18). 
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Let  us  then  take  the  subsystem  transforms 
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One  then  finds  the  following  transformed  problem 
(here  * - complex  conjugate) : 


x(k+l,t)  - A(l)x(k,l)  + B(i)u(k,i)  + w(k,i)  (4.10) 
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